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Abstract. We report on the 2005 AIM workshop "Generalized Kostka Polynomials", which gathered 20 
researchers in the active area of q, t-analogues of symmetric functions. Our goal is to present a typical 
use-case of the open source package MuPAD-Combinat in a research environment. 

1. Introduction 

The project MuPAD-Combinat (sec http://mupad-combinat.sf.net/ and |4]), born in spring 2001 under 
the leadership of F. Hivert and N. Thicry, is an open source package for making algebraic combinatorics 
using the computer algebra system MuPAD (see |www . mupad . de j and [5] for more details). The main goal of 
this package is to bring an open source flexible and easily extensible toolbox for checking conjectures in a 
short programming time. This package contains a large collection of tools as implementations of classical 
combinatorial objects (partitions, Young tableaux, trees, . . .), computations in combinatorial Hopf algebras 
(in particular symmetric functions), manipulations of graphs, automata, . . . 

The study of symmetric functions is a major historical field in algebraic combinatorics ([5]) and some 
people are mainly interested in generalizations of Kostka polynomials. In section 2, we give briefly mathe- 
matical definitions about different objects we work on. We illustrate how a computer algebra system which 
contains classical combinatorial objects, evolutive implementation of symmetric functions and easy way to 
incorporate program written in C++ can be used for our research. 

In section 3, we explain the design of our implementation of symmetric functions and some of our technical 
choices. We also give some examples of computations and show some advanced functionalities as adding new 
bases on the fly (using different characterizations) or implementing new operators on symmetric functions. 
Finally, in section 4 we describe how we incorporate programs into a coherent design using MAPITL library. 
The technical concepts used in section 3 and 4 are essentially classical, but we want to stress on their 
integration into the package in order to have a powerfull tool for making efficient. The difficulty is to 
find the appropriate combination that yields an intuitive yet flexible and powerful research tool. We show 
that our choices are promising in a real situation of collaboration at a workshop on generalized Kostka 
polynomials in Palo Alto, California organized by the American Institute of Mathematics in July 2005 (see 
|http : //www . aimath . o rg/WWN/kostka/ for more informations on this event). 

2. Symmetric functions and Kostka polynomials 

2.0.1. Basic definitions. A symmetric polynomial in variables X — {x\, . . . ,x n } is a polynomial in X in- 
variant under permutations of variables. When is infinite, we call symmetric function such a polynomial. 
The set of symmetric functions with coefficients in C(t), denoted A t , is a graded algebra with respect to the 
degree of polynomials, i.e 

A t = ®„> A?. 

For all n > 0, the dimension of A" is the number of partitions of n (a partition of a positive integer n, 
written A h n, is a decreasing sequence of positive integers with sum n). One main basis of this algebra is 
constituted with monomial functions defined for all partitions by A = (Ai, . . . , A ra ), by 

veo(x) 
l 
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where 0(A) represents the set of all the permutations of A. Another interesting basis consists of the symmetric 
powersums, defined for all partitions A by 

n 

P\(X) = pai • • -Px n where Pa« ='Y] x j i - 

3=1 

A scalar product on A t can be uniquely defined by 

(px,P»)=S^l[(m t )U m ^\ 

where mj(A) represents the multiplicity of part i in partition A. Applying Gram-Schmidt process of orthonor- 
malization on the monomial basis with respect to this scalar product yields us the basis of Schur functions 
(sa)a- These functions are intensively studied (from an algebraic and combinatorial point of view) and one 
of the beautiful results is the Littlewood-Richardson rule, a combinatorial interpretation of the product of 
two Schur functions. 

2.0.2. Kostka polynomials. One can introduce a t-deformation of the previous scalar product by 

IW 1 

(PA,PA = 5 x ^l[(m % y. J] I — 

i>l i=l 

The orthogonalization of the monomial basis with respect to this scalar product defines the Hall-Littlewood 
functions P\(X;t). There exist two other families of Hall-Littlewood functions: on the one hand, the 
Q\(X; t) which are the dual elements of P\(X; t) with respect to the scalar product ( , ) t and on the other- 
hand Q x (X;t) which are the dual elements of P\{X;t) with respect to (, ). The expansion of the Q x on 
Schur functions is an algebraic way to define the Kostka polynomials K\ tfi (t) 

Q' x (X;t)= J2 KxAt)^- 

These polynomials have also a combinatorial interpretation, using the charge on semi-standard Young 
tableaux in [7], or the rigged configurations introduced by Kerov, Kirillov and Reshetikhin in [5]. 

2.0. 3. Generalizations of Kostka polynomials. Using fc-ribbon tableaux introduced by Lascoux, Leclerc and 
Thibon in [B], we define for each positive integer k and each partition a particular symmetric functions 
H x k \X;t). Their expansion on Schur functions gives us a way to define an increasing filtration of Kostka 
polynomials K^' (t) 

fff (X;t)= £4>)V 

M h|A| 

The expansion of K x k ^{t) on the monomial basis is a way to define unrestricted generalized Kostka polyno- 
mials. In order to generalize Kostka polynomials, there exist other combinatorial ways (with unrestricted 
rigged configurations recently introduced by L. Deka and A. Schilling in pQ) and algebraic ways (using the 
theory of crystal bases for quantum groups of type A n ) to generalize Kostka polynomials. Other generaliza- 
tions are given by M. Zabrocki using creation operators in [16 and by L. Lapointe and J. Morse introducing 
a i-deformation of fc-Schur functions in [S]. An interesting problem, explained in [10] . is to show that these 
generalizations coincide in some particular cases. 

3. Implementation of symmetric functions 

3.1. Design goals. Symmetric functions can be represented in many different ways, and in particular in 
different basis (powersum, elementary, Schur, monomials,. . .). As usual in computer science, it is essential at 
each point to use the appropriate representation, both for efficiency and interpretation of the results. What 
make the situation specific is the number of those representations. In particular, it is neither practical nor 
sometimes possible to implement explicitely all conversions. Instead we want to be able to only implement 
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a few and deduce the others by compositions or linear tranformations such as inversion or transposition. 
Furthermore to thratend, the user and also programmers should not need to know which conversion are 
implemented because this information is too volatile. For example, the addition of symmetric function not 
given in the same basis is possible 

» S::s([2,l]) + S : : QP( [2, 1] ) + S::p([2,l]); 

(t + 4) m[l, 1, 1] + (t + 1) m[3] + (t + 3) m[2, 1] 

The monomial basis has been choosen by the system because it minimizes the number conversions (the cost 
of the conversion is not taken into account). The same process accurs for making conversion between two 
bases as the expansion of the Hall-Littlewood function Q 311 on the monomial basis 

» S: :m(S: : qP ( [2 , 1] ) ) ; 

(t + 2) m[l, 1, 1] + (t + 1) m[2, 1] + t m[3] 

We also consider operators on symmetric functions. In most cases, they are easier to define (and consequently 
to implement) on one of the bases but not on all of them. Consequently, the system uses implicit conversions 
between bases in order to apply operators on any given basis. The main technologies used are linear algebra 
and overloading mechanisms which are not essentially new, but a great effort is made in order to make 
manipulation intuitive. 

3.2. General overview of the implementation. In our design, for each basis of symmetric functions there 
is a domain (in the category Cat : : GradedHopf AlgebraWithBasis) which represents the space of symmetric 
functions expanded on this basis. For example, here is the example of implementation of the complete basis 

domain SymT: : complete (R: D0M_D0MAIN) 
inherits SymT :: common (R) ; 

category Cat: : GradedHopf AlgebraWithBasis (R) , Cat: :CommutativeRing; 
info_str := "Domain for symmetric functions expanded on complete basis"; 

basisName := hold(h) ; 

// Implementation of multiplication (complete basis is a multiplicative basis) 
mult2Basis := dom: : term@revert@sort@_concat ; 

end_domain: 

The domain Sym of symmetric functions, representing symmetric functions in whatever representation, is 
in category Cat: :Hopf AlgebraWith SeveralBases. This category helps us in defining implicit conversions 
between bases 

domain Sym(R=Dom: :ExpressionFieldO) 
inherits Dom: :BaseDomain; 

category Cat: :Hopf AlgebraWithSeveralBases CR) , Cat: :CommutativeRing; 

// Each domain corresponding to a basis is declared 
h := SymT: : complete Cdom: :coeff Ring, Options); 

// Implementation of different conversions between bases (stored in a table) 

basisChangesBasis := 
table ( 

// Explicit combinatorial conversions 

(dom::QP, dom::s) = (part -> (_plus (combinat :: tableaux : :kostkaPol (mu, part, dom::vHL) 

* dom : : s (mu) $ mu in 

combinat : : partitions : : list (_plus (op (part) ))))), 

(dom::McdP, dom::m) = (part -> ( (dom : : GramSchmidt (dom : :m, _plus (op(part) ) , 

dom: : scalartq) ) [op(part)] ) ) , 

// Dual conversions and inverse conversions 

(dom::s, dom::QP) = dom: : invertBasisChangeCdom: :QP, dom::s), 

(dom::s, dom::m) = dom: :transposeBasisChange(dom: :h, dom::s, dom::s, dom::m), 

) 

Note that only some conversions are implemented in this table; the overloading mechanism is in charge of 
finding the shortest number of intermediate conversions needed (it doesn't take into account the cost of each 
conversion). In order to use (q, i)-deformation of symmetric function, we can declare 
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» S :=examples : :SymmetricFunctions(Dom: :ExpressionFieldWithDegreeOneElements C [t ,q] ) , 

vHL=t , vMcd=q) ; 

3.2.1. Adding new bases on the fly. In order to add new basis on the fly, we define a generic domain. 

domain SymT: :NewBasis(R: D0M_D0MAIN, DomName : DDM_STRING) 
inherits SymT :: common (R) ; 

category Cat: : GradedHopf AlgebraWithBasis CR) , Cat: :CommutativeRing; 
info_str := "Domain for symmetric functions expanded on a new added basis"; 

basisName := text2expr (DomName) ; 

end_domain: 

In order to add a new basis named E for example, we can use 

» B := S: :newBasis(S: :coeffRing, "E"); 

and we define the change of bases we want. Let suppose that the change of basis between B and monomial 
basis corresponds to the function testChange 

>> S : : declareBasisChangeBasis C (B , dom: :m, testChange); 

>> S :: declareBasisChangeBasis Cdom : :m, B, dom: : invertBasisChange CB , dom::m)); 

3.2.2. Adding new operators. If we want to define new operators on symmetric functions, we only have 
to define their action on a particular basis, of course for efficiency. First, we declare our operators as an 
overloaded operator 

» newOp := operators: : overloaded C 

(x,y)->error( "Don't know how to compute the newOp on " . expr2text Cdomtype Cx) ) ) , 
Name="newOp") ; 

and by assuming that the action on the Schur basis is given by a function f , we declare 

» operators: : overloaded: : declareSignature C 
S: :newOp, [S : : s , D0M_INT] , 
S : : s : :moduleMorphism(f , dom: : s) ) ; 

and you can apply this operator on any basis due to the overloading mechanism. 

4. Integration of others programs written in C++ 

In July 2005, the American Institute of Mathematics organized in Palo Alto, California, a workshop on 
the generalized Kostka polynomials under the leadership of A. Schilling and M. Vazirani. During problem 
sessions, MuPAD-Combinat was put to use for testing conjectures. As usual in algebraic combinatorics, 
computations required the combination of preexisting combinatorial and algebraic functionnalities (as pro- 
vided by MuPAD-Combinat) with new combinatorial features (namely a highly technical bijection between 
fc-tuples of Young Tableaux and unrestricted rigged configurations) . Thanks to a dynamic module we could 
reuse a pre-existing robust C++ implementation of this bijection written by L. Deka; this allowed us to start 
manipulating large examples quickly. In general, dynamic modules permit us to reuse C++ code with two 
goals in mind: to avoid reimplcmcnting nontrivial and tested code, and to get quicker computations than in 
pure MuPAD language. In section 4.1, we describe the implementation of combinatorial objects in MuPAD, 
and in section 4.2 we present our technical choices for a seamless integration of a combinatorial bijection 
implemented in C++. 

4.1. Implementation of combinatorial objects in MuPAD-Combinat. We implement each combina- 
torial class as a domain in the MuPAD category Cat: : CombinatorialClassWith2DBoxedRepresentation. 
This category provides, among other things, a pretty printing method using ASCII characters. Let us 
illustrate this design in the case of skew riggings implemented in the domain combinat : : skewRiggings. 

domain combinat : : skewRiggings 
inherits Dom: :BaseDomain; 

category Cat: :CombinatorialClassWith2DBoxedRepresentation; 
axiom Ax : : canonicalRep ; 

info_str := "Combinatorial class for rigged skew partitions"; 
end_domain : 

We implement the constructor combinat: : skewRiggings : :new which builds an object from the list of its 
operands: 
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» a:=combinat: : skewRiggings ([ [] , [[[0], [0, 1]], [[""], [0, 1]]]]); 

H h 

10 10 1 
H + 

I 10 1 

H + 

The first element of the internal representation is the type 

» a[0] ; 

combinat : : skewRiggings 

Ribbon rigged configurations are particular sequences of skew riggings, implemented as plain lists of typed 
MuPAD objects. We implement them in the following domain 

domain combinat: : riggedConf igurations : :RcRibbonsTableaux 
inherits Dom: :BaseDomain; 

category Cat: :CombinatorialClassWith2DBoxedRepresentation, 

// Elements of this domain are represented using a plain MuPAD data structure 
Cat : :FacadeDomain(D0M_LIST) ; 

info_str := "Combinatorial class for rigged configurations"; 

end_domain : 

Here is an example of the bijection applied on the set of all 3-ribbon tableaux of shape (432) and evaluation 
(111) 

>> rc := map(combinat : : ribbonsTableaux : :list( [4,3,2] , [1,1,1] ,3) , 

combinat : : riggedConf igurations : :RcRibbonsTableaux : : f romRibbonTableau) ; 

II III 

| | + + + + + + + + | 

I I 1110 1,|0|0|0 0,1 I ||, 

I — + + + + + + + + — 



I III 

I + + + + + + + + I 

I 1010 1,101010 0,1 I I I, 



I +— + II II 

I +— + 10 10 +— +— + I I 

I 10 10 , +— + ,1 I III 

I +— + I 10 1 +— +— + I I 



» a:= rc[l] ; 

+ + 

I III 

| + + + + + + + + | 

I 1110 1,101010 0,1 I I I 

— + + + + + + + + — 

» op (a) [2] [0] ; 

combinat : : skewRiggings 



L. Deka and A. Schilling introduced in QQ a new kind of rigged configurations, namely the unrestricted 
ones. They were implemented as an independent CH — h program (file: FromOneCrystalPath.ee, headers: 
FromOneCrystalPath .h). On each rigged configurations (ribbons one and unrestricted one) we can compute 
a statistic. The interesting question is to find, in a special case, a bijection which preserves the statistic 
between these two kinds of rigged configurations. In order to manipulate these two objects in a single 
program, we decided in collaboration with A. Schilling and L. Dcka, to also integrate this C++ program into 
MuPAD- Combinat. In order to make the integration of this version of rigged configurations in a transparent 
way for the user, we kept the same design we used for ribbon rigged configurations. 
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4.2. Implementation of combinatorial objects in a dynamic module using the MAPITL library. 

We describe now the integration of the previous C++ program using MAPITL library (MuPAD Application 
Programming Interface Template Library (included in MuPAD-Combinat). This library provides 

• Wrappers to use MuPAD lists as standard containers 

• Easy C++ < — > MuPAD conversions with one single overloaded template for each directions: 

- C++ object — > MuPAD object: CtoM(c) 

- MuPAD object — > C++ object: MtoC(c) 

This includes conversions to/from containers and recursive containers. 

The C++ program computes the rigged configurations corresponding to a list of Young tableaux also 
called path which is denoted by the class path_class . We want to call from MuPAD 
void path_class: : buildj:igged_f or_path using the procedure 
combinat: : RiggedConf igurations : :newRiggedConf igurations. 
The building of the interaction is divided into three steps which are realized in the file 
RiggedConf igurationsPaths .mcc 

• convert a list of Young Tableaux in MuPAD to a C++ object of the class path_class 

• call the function void path_class : :buildjriggedJor_path 

• convert the result into a MuPAD object of type D0M_LIST 

Original files arc FromOneCrystalPath. cc containing declarations of different classes. We now explain some 
parts of the file RiggedConf igurationsPaths .mcc which is entirely given in appendix. The beginning of 
the file is the inclusion of the header file of the independent program and MAPITL 

#include<vector> 

#include "FromOneCrystalPath . h" 

#include " MAPITL. h" 

#include "MAPITL_tmpl . h" 

4.2.1. Conversion MuPAD to C+ + . 
MFUNCC newRiggedConf igurations , MCnop ) 
{ 

MFnargsCheckO) ; 
MFargCheck ( 1 , D0M_INT) ; 
MFargCheck(2, D0M_INT) ; 
MFargCheckO, D0M_LIST) ; 

n = MtoC<int>(MFarg(D) ; 

int path_len = MtoC<int> (MFarg(2) ) ; 

Cell input (MFargO) ) ; 

Cell *toto = input . toArray<Cell> ; 

> 

MFargCheck is the type checking of MuPAD arguments and MtoC permits to convert MuPAD objects into 
C++ ones. After we initialize an object of class path_class with values contained in toto 

path_class* input_path = new path_class (path_len) 
fordong a=0; a < input . size () ; a++) 

{ 

} 

4.2.2. Using the original function in C++. We can call now the original function in order to compute the 
bijection 

input _path->build_rigged_for_path() 

4.2.3. Conversion C++ to MuPAD. Next we create a C++ vector which contained operands of the corre- 
sponding into MuPAD object 

for(i=0; i<n; i++) 
{ 

vector<int> yyy; 
j = 0; 

«hile(input_path->rigged[i] [0] [j] != UNUSED) 
{ 

yyy .push_back(input_path->rigged[i] [0] [j]); 
yyy .push_back(input_path->rigged[i] [2] [j]); 
yyy .push_backCinput_path->rigged[i] [1] [j]); 
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} 

res[i] = yyy; 

} 

and we return the MuPAD list 

MFreturn(CelKres)) ; 

4.3. Compilation of the dynamic module. The next step is the compilation of the mcc file giving us the 
dynamic module RiggedConf igurationsPaths .mdm. We have to load this module from MuPAD using the 
function 

combinat : :RiggedConf igurationsPaths := 
procO 

save RiggedConf igurationsPaths ; 
begin 

if module : :which("RiggedConf igurationsPaths" ) = FAIL then 

userinfod, "Dynamic module RiggedConf igurationsPaths not available"): 

RiggedConf igurationsPaths := FAIL; 
else 

if traperror(moduleC"RiggedConfigurationsPaths")) <> then 

warning ("Error loading the dynamic module RiggedConf igurationsPaths :") ; 
lasterrorO ; 
end_if ; 
end_if ; 

RiggedConf igurationsPaths ; 
end_procO : 

We create the domain 

combinat: : riggedConf igurations : :RcPathsEnergy and procedure 

combinat: : riggedConf igurations : :RcPathsEnergy : : f romOnePath computes the bijection by calling the 
function implemented in the dynamic module 

combinat: : RiggedConf igurations : :newRiggedConf igurations 

>> a:= [combinat : : tableaux C [ [3,3]] ) , combinat :: tableaux ([ [2 , 2] ]) , combinat :: tableaux ([ [1 , 1] ])] ; 

— + + + + + + + + + — 

I 13 13 1,12 12 1,11111 I 

— + + + + + + + + + — 

We compute the bijection 

>> rc := combinat :: riggedConf igurations : :RcPathsEnergy :: fromOnePathCa) ; 

+— +— + | 
,110 1 I 

+— +— + | 

» op(rc[l]) [0] 

combinat : : skewRiggings 

The object computed using the C++ program is interfacing in a transparent way for the user who can use 
other MuPAD functionalities on the previous result. 

4.3.1. Acknowledgment. The author wants to thank A. Schilling and L. Deka for their collaboration on rigged 
configurations, M. Zabroki for its code on creating operators, and F. Hivert and N. Thiery for their support 
in the project MuPAD-Combinat. 
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Appendix : RiggedConfigurationsPaths . mcc 

#include<vector> 

#include "FromOneCrystalPath . h" 

#include "MAPITL.h" 

#include "MAPITL_tmpl . h" 

using MAPITL: :MtoC; 

using MAPITL: :CtoM; 

using MAPITL: : Container; 

using MAPITL : : SimpleCell ; 

using MAPITL: :Cell; 

using MAPITL : : Object InCell ; 

using MAPITL: : ArraylnCell ; 

using namespace std; 

MFUNCC newRiggedConf igurations , MCnop ) 
{ 

MFnargsCheckO) ; 

MFargCheckG, D0M_INT) ; 

MFargCheck(2, D0M_INT) ; 

MFargCheckO, DQM.LIST) ; 

n = MtoC<int>(MFarg(D) ; 

int path_len = MtoC<int> (MFarg(2) ) ; 

Cell input (MFargO) ) ; 

int i, j, k, trap, path_index, tblu_index, col; 

tmp = UNUSED; 
path_ index = 0; 
tblu_ index = 0; 

i = 0; j = 0; k = 0; col = 0; 1 = 0; 
initialize_lambda() ; 

path_class* input_path = new path_class (path_len) ; 
reset_tableau() ; 

Cell *toto = input . toArray<Cell> C) ; 
vector< vector<int> > resCn+1); 

fordong a=0; a < input . size () ; a++) 
{ 

Cell* alpha= (toto [a] ) . toArray<Cell> () ; 
int s = Calpha[0] ) . size () ; 

tblu_class *my_tblu = new tblu_class C (toto [a] ) . size C) ,s); 
my_tblu->tblu_id = tblu_index; 
tblu_index += 1 ; 

for (long d=0; d<toto [a] . size () ; d++) 
{ 

int* beta= (alpha [d] ) . toArray<int> () ; 

fordong b=0; b < s ;b++) 

1 

my_tblu->tb [d] [b] = beta[b]; 
my_tblu->tab_lambda[beta[b]-l] = 

my_tblu->tab_lambda[beta[b]-l] + 1; 

} 

} 
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input _path->path [path_ index] = my_tblu ; 
reset_tableau () ; 
path_index += 1; 

> 

input _path->build_rigged_f or _path() ; 

for(i=0; i<n; i++) 

{ 

vector<int> yyy; 
j = 0; 

while(input_path->rigged[i] [0] [j] != UNUSED) 
{ 

yyy.push_back(input_patli->rigged[i] [0] [j]) 
yyy.push_back(input_path->rigged[i] [2] [j]) 
yyy.push_backCinput_path->rigged[i] [1] [j]) 

} 

res[i] = yyy; 

} 

input_path->calculate_cocharge C) ; 
vector<int> stat; 

stat .push_backCinput_path->cocharge) ; 
res [n] = stat; 

MFreturn(CelKres)) ; 
> MFEND 



